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Abstract 


The leader RNA sequence was determined for two pig coronaviruses, tranmissi- 
ble gastroenteritis virus (TGEV), and porcine respiratory coronavirus (PRCV). 
Primer extension, of a synthetic oligonucleotide complementary to the 5’ end of 
the nucleoprotein gene of TGEV was used to produce a single-stranded DNA 
copy of the leader RNA from the nucleoprotein mRNA species from TGEV and 
PRCV, the sequences of which were determined by Maxam and Gilbert cleavage. 
Northern blot analysis, using a synthetic oligonucleotide complementary to the 
leader RNA, showed that the leader RNA sequence was present on all of the 
subgenomic mRNA species. The porcine coronavirus leader RNA sequences 
were compared to each other and to published coronavirus leader RNA se- 


The nucleotide sequence data reported in this paper have been submitted to the EMBL/Genbank/ 
DDBJ nucleotide sequence databases and have been assigned the accession numbers X52157, X52668. 
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quences, Sequence homologies and secondary structure similarities were identi- 
fied that may play a role in the biological function of these RNA sequences. 


Introduction 


Transmissible gastroenteritis virus (TGEV) and porcine respiratory coronavirus 
(PRCV) belong to the family Coronaviridae, a large group of pleomorphic envel- 
oped viruses with a positive-stranded RNA genome. TGEV causes gastroenteritis 
in pigs, resulting in a high mortality in neonates (1). PRCV was isolated in several 
European countries between 1984 and 1986 (2-4), does not cause diarrhea, and 
has been shown to replicate in the respiratory tract with little or no clinical signs, 
but is very similar antigenically and serologically to TGEV (2,4). Virions from 
both viruses contain two envelope glycoproteins of relative molecular mass (M,) 
200,000 (spike) and M, 28,000—31,000 (membrane protein) and a phosphorylated 
nucleoprotein of M, 47,000. cDNA probes to the structural protein genes of 
TGEV hybridized to the appropriate mRNA species of PRCV, suggesting a high 
degree of homology at the RNA level (unpublished data). 

Coronavirus proteins are expressed from a ‘‘nested’’ set of subgenomic 
mRNAs with common 3’ termini but different 5’ extensions. The sequence of 
each mRNA that is translated to produce viral proteins appears to correspond to 
the 5'-terminal region that is absent on the preceding smaller mRNA species. It 
has been shown for the coronaviruses, mouse hepatitis virus (MHV) and infec- 
tious bronchitis virus (IBV), the subgenomic mRNA species possess short 
‘leader sequences”’ at their 5’ ends. These sequences are not transcribed as a 
contiguous mRNA species, but are derived from the 5’ end of the genomic RNA 
and are probably joined to the 5’ end of each mRNA by a process of discontinuous 
transcription (S—9). The leader sequence appears to be produced by a mechanism 
termed leader-primed transcription, in which the leader RNA is transcribed inde- 
pendently, dissociated from the template, and then binds to the template (nega- 
tive-sense strand) at specific transcriptional start sites (10,11). The mechanism 
appears to involve the recognition of consensus sequences identified on the geno- 
mic RNA at those points corresponding to the 5’ ends of the subgenomic mRNAs. 
These consensus sequences may act as a binding site for the RNA polymerase- 
leader complex (7-9, 12-14). It has been previously postulated that a heptameric 
sequence, ACTAAAC (15-17), or a hexameric sequence, CTAAAC (18-20), may 
be involved in the binding of the TGEV RNA polymerase leader. 

In this paper we describe the elucidation of the leader RNA sequences from 
the porcine coronaviruses TGEV and PRCV, the first leader sequence to be 
described from the TGEV serogroup of coronaviruses. Comparison of the leader 
RNAs of TGEV and PRCV with published leader RNAs of other coronaviruses 
was used to identify areas of conserved sequence and potential secondary struc- 
ture that may be involved in the transcription of coronavirus subgenomic mRNA 
species. 
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Materials and Methods 
Preparation of viral RNA 


Confluent cultures of a pig kidney cell line LLC-PK1 were infected with a virulent 
British field isolate of TGEV strain FS772/70 or a British isolate of PRCV strain 
86/137004 at a MOI of 1-10 PFU per cell. After 2 hr at 37°C, the inoculum 
was removed and replaced with medium containing | pg/ml actinomycin D to in- 
hibit host-cell RNA synthesis (21). After a further 2-hr incubation, 25 pCi of 
[5,6-*H]uridine (Amersham International plc, TRK.410, 35-50 Ci/mM) was added 
per culture bottle and the cells were incubated for a further 5 hr. The cells were 
lysed with guanidinium thiocyanate, the RNA pelleted through 5.7 M cesium 
chloride and poly(A)-containing RNA isolated by poly(U) Sepharose affinity chro- 
matography, as described previously (21). 


Synthesis of oligonucleotide primers 


Two oligonucleotides were synthesized by the phosphoramidite method using an 
Applied Biosystem 381A synthesizer. One oligonucleotide, oligo 38 (5'’-TGGATT- 
CATCCCCCCAACTA-3’), was complementary to the nucleoprotein gene 22 bp 
downstream from the initiation ATG codon (15), as shown in Fig. 1, and was 
used for primer extension. The second oligonucleotide, oligo 58 (5'-AGAGATA- 
TAGCCACGCTACACTCACTTTAC-3’), was complementary to the 5’ end of 
the leader RNA (Fig. 1) and was used for Northern blot analysis of viral mRNA. 


Primer extension of nucleoprotein mRNA 


Gel-purified oligo 38 (500 ng) was 5’'-end-labeled (22) using 20 U of T, polynucleo- 
tide kinase (Gibco-BRL, Paisley) and 20 Ci [y-**7P] ATP (Amersham International 
plc, PB 10168, 3000 ci/mM. Poly(A)-containing RNA (1.5 jg) isolated from 
TGEV- and PRCV-infected cells was resuspended in water and heated at 60°C 
for 3 min. A further incubation was carried out using the two mRNA preparations 
in 27 wl reaction volumes containing 40 U of RNasin (Promega Biotec, Liver- 
pool), 50 mM Tris-HCl (pH 8.3), 10 mM MgCl, 35 mM KCI, 30 mM 2-mercapto- 
ethanol, 3 mM dithiothreitol, 4 mM dNTPs, 5’-end-labeled oligo 38 (120 ng), and 
21 U of AMV reverse transcriptase (Super-RT, Anglian Biotech Ltd, Colchester) 
for 90 min at 42°C. Formamide dye (80% formamide, 10 mM NaOH, 1 mM 
EDTA, 0.1% xylene cylanol blue, 0.1% bromophenol blue) was added and the 
mixture boiled for 3 min and electrophoresed on a 40 cm buffer gradient sequenc- 
ing gel (23). The wet gel was autoradiographed for 1 hr to locate the primer- 
extended products, which were excised from the gel. The labeled fragments were 


292 PAGE, BRITTON AND BOURSNELL 


ol igo 58 


L CCTTTTAAAGTAAAGTGAGTGTAGCGTGGCTATATCTCTTCTTTTACTTTAACTAGCTTT 
L GTGCTAGATTTTGTCTTCGGACACCAACTCGAACTAAACTTCTAAATGGCCAACCAGGGA 
G AGTGAGCAAGAAAAATTATTACATATGGTATAACTAAACTTCTAAATGGCCAACCAGGGA 
be 
S E Q@ E K L LH M V * MAN Q G 


‘ Oligo 38 


G CAACGTGTTAGTTGGGGGGATGAATCCACCAAAATACGTGGTCGCTCCAATTCC 


Q R V S W GOD ES T K I R GR SS WN S 


Fig. 1. Alignment of sequences from the 5’ end of TGEV nucleoprotein mRNA (L) and the corre- 
sponding region of the TGEV genome (G). The positions of homology between the genome sequence 
and the leader sequence are shown; the point of divergence between the sequences is identified by 
an arrow. The positions of the synthetic oligonucleotides oligo 38, used for primer extension, and 
oligo 58, used to probe viral RNA for the presence of the leader RNA sequence, are shown. The 
direction of the arrows indicate that both of the oligonucleotides were the complement of the se- 
quences shown. The amino acid sequences below the genomic sequence represent the carboxyl 
terminus of the membrane protein (16) and the amino terminus of the nucleoprotein (15). The postu- 
lated TGEV RNA polymerase-leader complex binding site or consensus sequence is underlined. 


eluted from the polyacrylamide gel and chemically cleaved (24). Samples of the 
cleaved products from each of the primer extended products were electropho- 
resed on 6% polyacrylamide gels at 35 W constant power for two different lengths 
of time. 


Northern Blot Analysis 


TGEV and PRCV poly(A)-containing RNA was glyoxylated and separated on a 
1% agarose gel (22). The RNA was transferred onto Biodyne A membranes (Pall 
P/N BNNG3R 1.2 wm, Gallenkamp) in X20 SSC (X1 SSC = 0.15 M NaCl, 0.015 
M trisodium citrate, pH 7.0) for 18 hr and baked at 80°C for 2 hr. The membrane 
was boiled in 50 mM Tris-HCl pH 8.0 for 5 min to remove glyoxal groups from 
the RNA and prehybridized in the presence of 50% formamide for 6 hr at 42°C 
(15). The viral mRNA species were hydribidized with **P-labeled oligo 58 in the 
presence of 50% formamide for 18 hr at 42°C. The membrane was washed four 
times in X2 SSC containing 0.1% NaDodSO, for 15 min at room temperature and 
autoradiographed. 
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Results 
Sequence of the nucleoprotein mRNA leader 


Following primer extension, using oligo 38 at the 5’ end of the nucleoprotein 
gene from the porcine coronaviruses TGEV and PRCV, labelled fragments of 
approximately 140 bases were produced and purified from gels. Larger molecular 
weight species were also observed (data not shown) in minor amounts, presum- 
ably corresponding to read-through sequences upstream of the nucleoprotein gene 
primed from the larger mRNA species. The nucleotide sequences of the two 
fragments, determined by chemical cleavage, were identical. The resulting nucle- 
otide sequence of the TGEV leader RNA sequence is shown in relation to the 
TGEV nucleoprotein gene in Fig. 1. The leader RNA sequence diverges from the 
genomic sequence 15 bp upstream of the nucleoprotein gene, corresponding to 
the first nucleotide of the membrane protein gene stop codon (16), indicating a 
length of 91 nucleotides of unique sequence (Fig. 1). The 91 nucleotide leader 
sequence of TGEV and PRCV has a low content of G (18%) and C (20%), and a 
high A (22%) and T (40%) content, with 20% of the T residues grouped in three- 
to four-nucleotide motifs (Fig. 1). These values are similar to those observed from 
the TGEV genome so far sequenced, except that the values for A (30.5%) and T 
(32.1%) are more similar on the genome than on the leader sequence. 

Analysis of the TGEV nucleoprotein nucleotide sequence (15) revealed a poten- 
tial RNA polymerase-leader complex binding site. The site, ACTAAAC, is seven 
nucleotides upstream of the nucleoprotein initiation codon and has also been 
found to precede all the TGEV structural protein genes and two of the three 
potential genes shown to be at the 5’ end of mRNA species (15-17). This consen- 
sus sequence is found two nucleotides downstream of the nucleotide where the 
leader RNA and TGEV genomic sequences diverge, indicating that this sequence 
is involved in the leader-primed transcription of TGEV mRNA molecules. As can 
be seen from Fig. 2, 4 of the 6 mRNA species from the FS772/70 strain of TGEV 
have the sequence AACTAAAC, of which the 5'-end adenosine residue is the 
next base down from the divergence point. In fact, the consensus sequence at 
the spike/ORF1-ORF2 gene junction has the sequence GAACTAAAC and at the 
NUC/ORF4 gene junction has the sequence CGAACTAAAC, indicating that the 
region of the leader sequence 5’ to the homology motif, ACTAAAC, may vary 
between 89 and 91 nucleotides depending on the TGEV gene. 

Computer analysis has also detected a homology between the leader RNA 
sequence and the 5’ end of the negative strand (i.e., the reverse complement of 
the noncoding region at the 3’ end of the positive strand). This is shown in Fig. 
3. The nucleotides on the leader RNA sequence, bases 84-99, and on the negative 
strand, bases 136 to 152 counting from the first base after the poly(A) tail, have 
an overall homology of 82% and include the sequence CTAAAC, which is part 
of the postulated TGEV RNA polymerase-leader complex binding site. This is 
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Gene Junction Intergenic Sequences 


POL/SPIKE TAAGTTACTAAACTTTGGTAACCACTTCGTTAACACACCATG 


Spike/ORF1-ORF2 TTAAGAACTAAACTTTCAAGTCATTACAGGTCCTGTATG 


ORF2/ORF3 GGCGGTTCTAAACGAAATTGACTT AAAAGAAGAAGAGGGAGACCGTACCTATG 
ORF3/MEM GTTTGAACTAAACAAAATG 
MEM/NUC GGTATAACTAAACTTCTAAATG 
NUC/ORF4 TAACGAACTAAACGAGATG 
LEADER ACTCGAACTAAAC 
87 99 


Fig. 2. Comparison of the TGEV (strain FS772/70) gene junctions with the sequences immediately 
5' of the consensus sequence and part of the TGEV leader RNA sequence. The positions of the 
ACTAAAC consensus sequences and any identical bases 5’ to this sequence present on the leader 
RNA sequence are double underlined. The initiation codon of the gene immediately downstream of 
the consensus sequence is underlined. The sequences of the spike/ORF1-ORF2 and ORF2/ORF3 
junctions are taken from (17); the ORF3/MEM junction is taken from (16); and the MEM/NUC and 
NUC/ORF4 junctions are taken from (15). The POL/SPIKE junction sequence is from unpublished 
work. POL = polymerase; MEM = M, 29 459 glycoprotein; NUC = nucleoprotein. 


very similar to the observation for IBV (25) involving sequences present at the 
5’ end of the IBV genome, and on the IBV leader RNA sequences, with the S’ 
end of the IBV negative strand. The homology observed included the sequence 
CTTAAC, which is part of the postulated IBV RNA polymerase-leader complex 
binding site CT(T/G)AACAA. 


Northern blot analysis of TGEV and PRCV mRNA subgenomic species 


An oligonucleotide, oligo 58, was synthesised that was complementary to the 5’ 
end of the TGEV and PRCV leader RNA sequences (Fig. 1). The oligonucleotide 


TGEV GENOME (-) 115 BASTTACT ARDS ST CAAGUD NES? COAAATICARE PET BARC 152 


wee eee 
eeeses 


TGEV LEADER 62 TGCTAGATTTT. GTCTTCGCACA. CCAACTCCAA. CTAAAC 99 


Fig. 3. Comparison between the reverse complement of the 3’ end of the TGEV genome (i.e., the 5’ 
end of the negative strand (15) and part of the TGEV leader RNA sequence. Colons show identical 
bases. The single dots in the sequences are padding characters inserted to achieve optimal alignment. 
Part of the postulated TGEV RNA polymerase-leader complex binding site sequence or consensus 
sequence is double underlined. 
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mRNA 
species 


ani, G 


— ~ —__——§ 3.65 
, =: 
eo =a — 26 


Track 1 2 38 


Fig. 4. Autoradiograph of the mRNA species from two strains of PRCV and TGEV northern blotted 
and probed with the 5’ end-labeled synthetic oligonucleotide, oligo 58 (see Fig. 1). Track 1 PRCV is 
strain 86/135308, track 2 is TGEV strain FS772/70, and track 3 is PRCV strain 86/137004. The ORF1/ 
ORF2 mRNA species from PRCV migrates faster than the corresponding TGEV species, relating to 
a size of 3.65 kb (34). 


was end-labeled and used to probe TGEV and PRCV mRNA species that were 
Northern blotted onto Biodyne membranes. As can be seen from Fig. 4, the 
labeled probe hybridized to all of the TGEV and PRCV mRNA species. The 
intensity of the bands corresponding to labeled probe hybridized the spike mRNA 
species, and genomic RNA was lower than that observed for the smaller mRNA 
species due to less of these larger species being isolated from the poly(U) Sepha- 
rose column used in the isolation of mRNA. The fact that the probe hybridized 
to all of the mRNA species showed that the leader RNA sequence was present 
on the other RNA molecules of TGEV and both strains of PRCV was not unique 
to the nucleoprotein mRNA species. 
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Homology to other coronavirus leader RNA sequences 


The two porcine coronavirus leader sequences were identical, indicating that the 
two viruses probably use the same RNA polymerase-leader complex binding site, 
ACTAAAC, for the synthesis of subgenomic mRNA species. The SEQHP com- 
parison program of the Los Alamos (26) package was used to compare the leader 
RNA sequences determined in this paper and those published for five other coro- 
naviruses belonging to two different serogroups. The sequences were compared 
from the 5’ ends to the point of divergence from the genomic sequences. The 
percentage homologies, Table 1, were expressed as the number of bases matched 
to the longer of the two sequences being compared. The homology of the leader 
sequences fell into three groups. Leader RNAs from coronaviruses belonging to 
different serological groups had homologies in the region of 35-40%. Serologically 
related viruses like human coronavirus (HCV) (strain OC43) and MHV (strains 
AS9 and JHM) have about 60% homology. The third group involved different 
strains of MHV, AS59, and JHM, which showed a homology of 91%. This observa- 
tion indicates that TGEV and PRCV, which have a homology of 100%, are proba- 
bly different strains of the same virus or that PRCV has very recently diverged 
from TGEV. 

In order to identify common areas of homology, the leader RNA sequences 
from seven coronaviruses were aligned. As can be seen from Fig. 5, these fell 
into two groups. One group consists of MHV (strains A59 and JHM) with HCV 
(OC43), which have a fairly high degree of homology along their lengths. The 
other group consists of TGEV and PRCV (not shown on the diagram) with HCV 
(229E) and IBV, which have high homologies at their 3’ ends and areas of homol- 
ogy at their 5’ ends. There are good homologies towards the 3’ ends, involving 
the postulated RNA polymerase-leader complex binding sites and sequences up- 
stream of these sites, between the groups, but very little if any homology between 
the 5’ ends. 


Table 1. Comparison of coronavirus leader RNA sequences 


Percentage homology of leader RNA sequences* 
PRCV HCV HCV MHV MHV IBV 
(86/137004) (229E) (OC43) (A59) (JHM) (Beaudette) 


TGEV (FS772/70) 100 36 39 35 39 40 
HCV (229E) 100 45 42 41 52 
HCV (0C43) 100 62 67 38 
MHV (A59) 100 91 49 
MHV (JHM) 100 44 
IBV (Beaudette) 100 


4The percentage homologies were determined from the 5’ end of the leader RNA sequences to the 
points of divergence from the genomic RNA sequence. All values are calculated using the longest 
sequence of the pairs. 
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229E  CTTAAGT......... ACCTTATCTATCTACAAATAG. AAAAGTTGC. . .TTT.T. .TAGACTTTGTGTCTACTT. ...-...+20.040- CTAAAC 
TGEV act. GAGTGTAG. . - CGTGGCTATATCTCTTCTTTTACTITAACTAGC _. sTTTGIGCTAGATITICT.... CTTCGGACACCAACTOGAACTAAAC 
IBV ACTTAAG. ATAGATATTAATATATATCT aoe ATTAC. - -ACTAGC. s Ser TCOECTAGATERTTIK Mg SERTS secs As AAC. AAAA 
0ca3 TTGTGAGCGAAGTTC. catecetacar. cooceTTccActcat ox! ctorrertacatetrttt ....CTAAT....CTAA.T... .CTAAAT 
A59 TATAAGAGTGATTGGCGTCCCTACGTACCCTCT . CAACTCTAAAACTCTTG. TAG. TITAAAT.. ..CTAAT. << .ceceeeeees CTAAAC 
JH TATAAGAGTGATTGGOGTCCGTACGTACCCTCT .CTACTCTAAAACTCTTG. TAG. TTTAAAT... .CTAAT. . . .CTAA.T. .. .CTAAAC 

kk kk * . 4 the at kk ReeKKKKK kk kkkKKK 

* * tt * ake « ke kk kK 


Fig. 5. Alignment of leader RNA sequences from nucleotide 14 in the TGEV sequence with those 
from other coronaviruses. PRCV is not included on the diagram, as it is identical to TGEV. Colons 
between pairs of sequences show identical bases. The top row of asterisks marks positions in which 
four or more bases are identical, and the bottom row of asterisks marks bases completely conserved 
between all seven sequences. The dots within the sequences are padding characters inserted to achieve 
optimal alignment. The sequences of the leader RNA sequences were taken from TGEV (FS772/70) 
and PRCV (86/137004), this paper; HCV, strain 229E (31) and strain OC43 (35); MHV, strain A59 (7) 
and strain JHM (13); avian, IBV strain Beaudette (9,25). 


Prediction of secondary structures 


As seen from Fig. 5 simple alignment did not reveal very much information about 
the homologies of the leader RNA sequences from the different coronaviruses, 
except at the 3’ ends involving the consensus sequences. In order to identify any 
potential similarities in these sequences, the secondary structure of the RNA 
sequences in Fig. 5 were analyzed. Potential secondary structures of the leader 
RNA sequences were determined using the computer program FOLD (27) from 
the UWGCG DNA analysis programs (28). The coordinates determined by the 
FOLD program were displayed graphically using the UWGCG program SQUIG- 
GLES. The potential secondary structures obtained were compared and, as can 
be seen from Fig. 6, the overall shape of these sequences are very similar, except 
for the avian coronavirus IBV. All the molecules appear to be composed of two 
stem-loop structures. The two MHV molecules are very similar in shape and, as 
seen from Fig. 5 and Table 1, are very homologous, 91%, at base sequence. The 
secondary structures of the coronavirus leader RNA sequences are probably 
influenced by their biological function, which results in the similarity of these 
potential structures. 


Discussion 


This paper presents evidence that the nucleoprotein mRNA species of TGEV 
and the closely related porcine respiratory variant of TGEV, PRCV, contain an 
identical leader RNA sequence of about 91 nucleotides. Sequencing studies on 
TGEV have shown that the heptameric sequence ACTAAAC occurs on the ge- 
nome upstream of the genes and is believed to be the binding site for the leader 
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PORCINE (TGEV) FS772/70 LEADER SEQUENCE HUMAN 229E LEADER SEQUENCE 
Length: 115 Energy: -22.4 Length: 80 Energy: -10.7 


20 


100 


80 
HUMAN OC43 LEADER SEQUENCE MURINE (MHV) A59 LEADER SEQUENCE 
Length: 91 Energy:-13.9 Length: 72 Energy: -10.1 


20 60 


40 


MURINE (MHV) JHM LEADER SEQUENCE CHICKEN (IBV) BEAUDETTE LEADER SEQUENCE 


Length: 77 Energy: -10.7 Length: 66 Energy: -12.1 
40 


60 
40 


Fig. 6. Comparison of the predicted secondary structures of the coronavirus leader RNA sequence 
as determined by the UWGCG programs FOLD and SQUIGGLES (28). The leader RNA sequences 
used were as described in Fig. 5. The numbers refer to the positions of the nucleotides. 


RNA primer. Northern blot analysis, Fig. 4, of TGEV strain FS772/70 and PRCV 
strains 86/137004 and 86/135308 RNA, probed with a complementary oligonucleo- 
tide to the leader RNA identified on the nucleoprotein mRNA species of these 
coronaviruses, showed that all of the subgenomic mRNA molecules contained 
the common leader RNA sequence. These observations support the finding ob- 
served for MHV that the synthesis of coronavirus mRNA species is primed from 
the negative-sense strand using a small RNA molecule derived from the 5’ end 
of the genomic RNA. This mechanism has been termed leader-primed transcrip- 
tion and involves not only the leader RNA primer, but also consensus sequences 
along the genome found upstream of the genes, which act as binding sites for the 
leader RNA primer. 

Comparison of TGEV and PRCV viral products has shown very little difference 
between the two coronaviruses, and until recently is was impossible to differenti- 
ate between the two viruses using antisera. PRCV is fully neutralized by antisera 
prepared against TGEV, and the majority of monoclonal antibodies (MAbs) raised 
against TGEV virion proteins cross-react with PRCV. However, MAbs, raised 
against antigenic determinants of the spike protein from either the virulent British 
isolate FS772/70 (29) or the avirulent Purdue strain of TGEV (30) have been 
identified that do not recognize PRCV. These observations and the fact that the 
leader RNA sequences from TGEV and PRCV are identical supports the evidence 
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that the two viruses are very similar and that PRCV may have evolved as a TGEV 
variant. 

Comparison of the TGEV leader RNA sequence with the genomic sequence 
upstream of the nucleoprotein indicates that the length of the unique sequence of 
the leader sequence is 91 nucleotides. The point of divergence is two bases up- 
stream of the ACTAAAC sequence, supporting the evidence that the TGEV RNA 
polymerase-leader complex binding site is ACTAAAC. Four out of the six mRNA 
species from the FS772/70 strain of TGEV have the sequence AACTAAAC, and 
the 5'-end adenosine residue is the next base down from the divergence point in 
the nucleoprotein mRNA (Fig. 2). The differences in the homologies between the 
leader RNA and sequences upstream of the consensus sequence on the genomic 
RNA may play a role in the levels of transcription of a particular mRNA species. 
The mRNA species of 3.0 kb has been shown to have an open reading frame at 
the 5’ end encoding a potential polypeptide of M, 9200 (17). This particular mRNA 
does not have the heptameric consensus sequence but has the hexameric 
CTAAAC sequence, and it is interesting to note that it is the least abundant 
TGEV mRNA species (observed from TGEV mRNA in total cell lysates). Hy- 
bridization of oligo 58 to the 3.0-kb mRNA species showed that this species does 
contain the TGEV leader RNA, confirming that it is a true mRNA species, even 
though it is the only TGEV species not to have the heptameric consensus se- 
quence. 

Comparison of the seven coronavirus leader RNA sequences against each other 
identified three groups (Table 1): non-serologically related viruses had about 
35-40% homology; serologically related viruses had about 60% homology; viral 
strains had about 90-100% homology. However, TGEV and HCV (229E) have 
been placed in the same serological group, but have only 36% homology within 
their leader RNA sequences, suggesting that the two viruses are not particularly 
related. TGEV and HCV (229E) have been shown to have 46% homology at the 
amino acid level within their derived nucleoprotein sequences (31), whereas the 
homology between the derived nucleoprotein amino acid sequences for different 
viruses within the MHV serological group are between 80% and 98% homology. 
This indicates that the serological grouping of coronaviruses is not a particularly 
useful test, as similar epitopes may exist on the viral structural proteins. Compari- 
sons of nucleic and amino acid sequences from the viruses will provide a more 
accurate method for grouping the viruses. It will be interesting to compare the 
leader sequences of bovine coronavirus (BCV), which is serologically related to 
HCV (OC43) and MHV (A59 and JHM), with feline infectious peritonitis virus 
(FIPV) and canine coronavirus (CCV), which are serologically related to TGEV, 
once their sequences have been determined. 

The large variation in sequence length and content made the alignment of the 
different leader sequences difficult. However, alignment of the six different coro- 
naviruses revealed that they fell into two groups. There appears to be some 
conservation of short sequence motifs between the seven leader sequences. To- 
ward the 3’ end of the sequences, a TAG motif is conserved in all the leaders, 
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followed by a string of Ts. In five out of seven of the sequences, this motif is 
TAGANNTT. About ten nucleotides downstream of this region is a conserved 
CT motif, which is followed by a series of nucleotides differing in number, de- 
pending on the coronavirus, followed by the postulated RNA polymerase-leader 
complex binding site. The largest number of nucleotides between the CT motif 
and the consensus sequence are found on TGEV and PRCV, the shortest is found 
on HCV (229E) and IBV. It is interesting to note that there is a five-base insert 
in MHV strain JHM when compared to MHV strain A59, which is also present 
in HCV (0C43) within this region. All the mammalian coronaviruses appear to 
have the motive CTAAAC, except HCV (OC43), which has CTAAAT. Recent 
sequence data suggest that coronaviruses FIPV and BCV have ACTAAAC as 
their mRNA consensus sequence. Upstream of the TAG motif there is an ACT 
motif occurring in six out of seven sequences. Toward the 5S’ end of the leader 
RNA sequences, the homologies are patchy and limited to short matches, occur- 
ring only between pairs of sequences. 

The area upstream of the consensus sequence has been suggested to be in- 
volved in the binding of nucleoprotein to the leader RNA sequence at nucleotides 
56-65 in MHV (32). It was suggested that mRNA species and genomic RNA form 
a complex with the nucleoprotein by the protein binding to or near the leader 
sequence attached to the RNA molecules (33). Secondary structure analysis of 
the leader RNA sequences showed that all the sequences except for IBV possess 
a putative double stem-loop structure (Fig. 6). In the case of the mammalian 
coronaviruses, the consensus sequences and upstream regions of homology are 
on the second stem-loop structure, leaving the possibility that the RNA-depen- 
dent RNA polymerase could interact with the first stem-loop structure. The IBV 
consensus sequence is present on the free 3’ end of the single stem-loop structure, 
possibly leaving the single stem-loop structure to interact with the polymerase. 
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